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We investigate the dynamics of entanglement between two continuous variable quantum systems. 
The model system consists of two atoms in a harmonic trap which are interacting by a simplified s- 
wave scattering. We show, that the dynamically created entanglement changes in a steplike manner. 
Moreover, we introduce local operators which allow us to violate a Bell-CHSH inequality adapted 
to the continuous variable case. The correlations show nonclassical behavior and almost reach the 
maximal quantum mechanical value. This is interesting since the states prepared by this interaction 
are very different from any EPR-like state. 
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I. INTRODUCTION 

Since the seminal paper Q of Albert Einstein, Boris 
Podolsky and Nathan Rosen (EPR) and the ingenius 
criticism by John Bell fi || we have begun to under- 
stand that quantum correlations can be stronger than 
ever describable by a classical theory based on local re- 
alism. No local realistic theory (LRT) can fully explain 
all statistical properties predicted by quantum mechan- 
ics for specific composite systems. The Hilbert space 
structure allows for states of quantum systems which we 
call, according to Erwin Schrodinger 14], entangled. It 
was proven that entangled spin-i particles show corre- 
lations that violate the so-called Bell inequalities |^|, |). 
They provide a boundary on correlations that still can be 
modelled within a LRT. In fact, for any composite sys- 
tem described by a pure entangled state in a finite Hilbert 
space a set of observables exists that leads to a violation 
of the Bell inequalities |jj . Hence there is a close relation 
between entanglement and violations of Bell inequalities 
pointing to the nonclassicality of specific quantum me- 
chanical states. 

It became, however, soon clearer that the situation 
is more complicated when it comes to mixed states de- 
scribed by density operators. Here we need a finer dis- 
tinction of entanglement. One finds states which are non- 
separable but still fulfill the Bell inequalities || . Certain 
classes of these non-separable states can still be sub- 
jected to "entanglement distillation" and the resulting 
state then leads to a violation ||. In particular, any 
non-separable state of two qubits can be distilled |l(J . 

But it was also realized jll| that a bipartite system 
can be in a bound entangled state which satisfies the 
known versions of Bell inequalities. Therefore we can 
say that Bell-type arguments can indicate nonclassical 
correlations of finite-dimensional systems but they can- 
not characterize all kinds of entanglement jjll . 

The situation is again different when we move to con- 
tinuous variable systems with infinitely large Hilbert 
space. Correlated continuous variable states can be used 
as a nonclassical resource for quantum computation and 
quantum communication |L3| [l4|, |l5) . In particular, it 
was possible to realize unconditional teleportation in the 



continuous case p6| . Hence also in this realm entan- 
glement acts as an essential ingredient for nonclassical 
tasks. Consequently, criteria have been proposed |l7) in 
order to decide whether a continuous variable state is en- 
tangled and protocols have been discussed which purify 
continuous entangled states |l^, [llj . 

Composite continuous systems can also violate Bell- 
type inequalities. One approach is to map the continuous 
spectrum on a discrete subset of possible measurement 
outcomes for which the usual form of Bell inequalities 
applies |2C], ^|. This direction has a very appeal- 
ing phase-space interpretation and it establishes links 
between nonlocality and the Schrodinger-cat gedanken 
experiment p3[ . However, it is also possible to see non- 
classical EPR-type correlations in quadrature-phase mea- 
surements [p4[ which avoid any discretization. 

In the present paper we shall study an elementary con- 
tinuous variable system: two trapped particles that can 
collide with each other. In particular we are interested in 
the dynamical creation of entanglement between the two 
sub-systems and we shall show that it leads to an almost 
maximal violation of a Bell-type inequality. Neverthe- 
less, the generated entangled state is completely different 
from an EPR state. It is actually a superposition of EPR 
states widely separated in phase space. 

The paper is organized as follows: In Sec. II we give 
a description of the physical model we investigate. The 
mathematical description in form of a Schrodinger equa- 
tion and its solution is given. In Sec. |n] we investigate 
the entanglement contained in the prepared states using 
an entropic measure. Finally, in Sec. [fy| we introduce 
dichotomic observables which allow us to violate a Bell- 
type inequality. Sec. [v| adds the final conclusions. 



II. MODEL 

A variety of quantum mechanical systems can be con- 
sidered which have continuous degrees of freedom. By 
this we mean that on such systems we can in principle 
measure certain observables with a continuous spectrum. 
Examples are a quantized mode of the electromagnetic 
field with its quadrature observables and the quantized 
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motion of massive particles with the corresponding posi- 
tion and momentum observables. 

Engineering entanglement between two such sys- 
tems and exploiting this entanglement for quantum in- 
formation purposes requires that the systems can be dis- 
tinguished via any degree of freedom which does not take 
part in the entanglement process. If we consider, for ex- 
ample, atoms and entangle their motion, they still can 
be distinguished via their internal electronic structure. 

With current experimental techniques it is possible to 
trap ultracold atoms in an almost perfect harmonic po- 
tential. Moreover, one can address single atoms by excit- 
ing their internal structure with the help of lasers. In 
general the atom-atom interaction is a highly compli- 
cated processs. We assume that for cold atoms with a 
large de-Broglie wavelength we can neglegt the details 
of the interaction potential and replace it by a pointlikc 
interaction. For a composite system consisting of only 
two particles with mass m the corresponding Hamilto- 
nian then reads 



H 



Y~ (Pi +Pl) + \ m u 2 (x\ + xl)+h~{8{xi-X2) (1) 



with trap frequency lj and interaction strength hj > 
for a pointlike repulsive potential. 

Certain aspects of the dynamics of this model have al- 
ready been studied p6[ [27[] . Here we shall concentrate on 
the dynamics of the entanglement induced by the delta 
interaction. In order to do this we first derive the ener- 
gies E and eigenfunctions |\&) of the time independent 
Schrodinger equation 



(2) 



in appropriate coordinates. If we use a dimensionless 
center-of-mass coordinate 



Vc = k c ■ - (xi + x 2 ) with k c = 
and a dimensionless relative coordinate 



y r = K r ■ [x\ — X2) with K r 



2mu) 



~2h 



(3) 



(4) 



as well as scaled energies E = E/huj and a dimensionless 
interaction constant 7 = "fK r /uj we get 



dyl + dy 2 r ^ c + Vr > 

+2(E-l5(y r ))]y(y c ,y r ) = 0. (5) 



This Schrodinger equation can be solved using a sep- 
aration ansatz. We present the calculations in Ap- 
pendix |X| We find the eigenvalues E c and eigenfunctions 
tl> c (E c ; y c ) for the center-of-mass motion as well as E r and 
tl) r (E r ;y r ) for the relative motion of the particles. The 
total energy then reads E = E c + E r . We can combine 



these solutions to get the full solution of the stationary 
Schrodinger equation, Eq. (H). 

The dynamics of the system for scaled times t = tot is 
then completely given by the wave packet 

^(y c ,yr,t) = ^2 a(E c ,E r )tp c (E c ;y c )^ r (E r ;y r ) 



E c ,E r 



x exp 



(-i(E e + E r )tj (6) 



with complex coefficients a(E Cl E r ) determined by the 
initial condition ty(y c ,y r ,t = 0). Rescaling and trans- 
forming back the coordinates yields 

X\ + x 2 



y r = K r (x\ -x 2 ),t = u)t) (7) 

for the particle coordinates x\ and x 2 at time t. 

In the following we shall concentrate on the entangle- 
ment between the two trapped atoms and in particular on 
its dynamical evolution due to the pointlike interaction. 



III. ENTANGLEMENT MEASURE AND 
INTERACTION 

Most examples of possible applications of entangle- 
ment in continuous variable quantum information is 
based on EPR-states |jj or in the optical case approx- 
imations of EPR-states using squeezed states mixed at 
a beam splitter jl^, |2^, The time evolution in our 
system now produces entangled states of the two parti- 
cles quite different from EPR-like states. We will see this 
by looking at the Wigner function of the relative motion. 
Nevertheless, within specific time intervals our particles 
enjoy a strong nonclassical correlation. In order to dis- 
cuss the strength of the particles' correlation and the non- 
local properties of the corresponding continuous- variable 
states we have to characterize this entanglement. In this 
section we will use a measure of entanglement based on 
the von-Neumann entropy, in Section IV we present a 



Bell-test of the nonlocality of the states prepared. 

A first hint of the atomic correlation produced by the 



Hamiltonian H , Eq. 
tanglement measure 
Neumann entropy 



can be given by a common en- 
. For pure states 1^) the von- 



S = -tr 



-red 



log Q 



rcd\ 



(8) 



gives a reasonable measure of entanglement when we in- 
sert the reduced density operator 



Q 



red 



tr|*)(*| 



(9) 



Here the trace is taken over one atomic subsystem. Hence 
in our case we obtain 



g Ted (x 1 ,x' 1 ,t)= dx 2 y(x 1 ,x 2 ,t)V*{x' 1 ,X2,t) (10) 
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for the state $>(x 1 ,x 2 ,t), Eq. (Q). 

In order to calculate the von-Neumann entropy of 
g rcd (xi : x[, t), Eq. (|8|), we have to diagonalize the re- 
duced density operator in order to find the eigenfunctions 
4>v(x\,t) which fulfill 



dx\ g rc (xi,x' 1 ,t)(j) v (x' 1 ,t) = \ u (t)<f> v (xx,t) (11) 



with eigenvalues A„(t). Note that these eigenfunctions 
form a complete and orthonormal set of functions for 
any time t. The integration, Eq. ([l0|), and the diagonal- 
ization, Eq. (|IT|), can be done numerically using discrete 
values for x%, X2 and x[. 

With help of the eigenvalues \ v (t) and eigenfunctions 
4> v {x\,t) we can write 

oo 

Q red (x 1 ,X , 1 ,t) = J2 X »(t)<t>»(xi,t) ( j>:(x' 1 ,t) (12) 
v=0 

and the von-Neumann entropy then reads 

oo 

S(t) = -J2Mt)^gX v (t). (13) 

i/=0 

Our question now is how this entropy of entanglement 
develops in time starting from an uncorrelated Gaussian 
product state 

/raw / muj , , ,t\ 

*( Xl ,x 2 ,t = o) = J— cxp (-—( Xl -R/2y) 

xexp(-^(x 2 + i?/2) 2 ) (14) 

for two particles with mean distance R. Hence the two 
atoms start symmetrically in the trap potential with 
no initial momentum and move towards the trap center 
where they collide. 

When we rewrite this state in the scaled center-of-mass 
coordinate, Eq. (||), and the relative coordinate, Eq. (||), 
it reads 



V M (y c ,y r ,t = 0) = -Lexpf-iy c 2 



, rxp | -I(y r _p) 2 ) . ( 1.51 



Here we have introduced the scaled distance parameter 
p = K r R which now completely characterizes the whole 
family of dimensionless initial wave packets, Eq. (15). 
Note further that these initial conditions mean that 
the center-of-mass motion parametrized by y c has been 
cooled to the ground state. Moreover, we recall by look- 
ing at Eq. (||) that the center-of-mass motion and the 
relative motion decouple. Hence, due to the chosen ini- 
tial conditions the complete dynamics takes place in the 



relative motion and our wave function has the product 
structure 

* W (lfci!fr,*) = -lj exp [~\vl\ ex P (-5*) 

x 4"Hy r ,t) (16) 

for all times t. Consequently, the characteristic features 
of the system are determined by the wave function ip r of 
the relative motion. 

First we will have a look at the short-time behavior of 
this wave packet. Fig. [I] shows the Wigner function 



W(y r ,p r ,t) = ^ /^expHK)^ 3 *^-^,*) 



(17) 



of the relative wavefunction ip r at times t = (a) and 
t = Air (b). Each time, the wave packet crosses the 6- 
peak, part of it is reflected coherently while the rest of 
it passes through the origin. This leads at time t — Air 
to the humps which are almost Gaussian and opposite to 
each other in the right part of Fig. [ll They clearly add 
up coherently as can be seen from the interference fringes 
in between. 

The state at time t = Air is certainly no longer separa- 
ble, but it is also much different from an EPR eigenstate. 
This can be seen if we look for the overlap with a prop- 
erly scaled EPR-state \Y r ,P c ) derived in App. ||. The 
main contribution to the overlap, Eq. (B6), 



(r r ,p c |*W(t)) 



,exp(-P c 2 ) (18) 



is given by the wavefunction ipl^ (Y r , t) of the relative 
coordinate. If we now look again at Fig. |l](b) we can 
see, that this overlap is not a sharp peak, but obviously 
delocalized. 

Let us understand further how the entanglement builds 
up in the state Eq. ([l6|) for short times. In fact we shall 
see that this is a steplike process: whenever the two par- 
ticles collide entanglement is added to the system. 

In order to see this we show in Fig. |^(a) the entangle- 
ment S(t), Eq. (||), as a function of scaled time t = uit 
together with the average interaction energy 



E int (t) 



S(xi -xi))= j(5(y r ) 



7 4 P \y r =0,t) 



(19) 



Here we have used Eqs. (Q) and (0). This energy will 
help us to identify the time intervals in which the inter- 
action takes place. We have plotted these quantities in 
Fig. |2| for time t ranging from to 2n. This is exactly 
the time for a complete round-trip in the case of no in- 
teration. 
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FIG. 1: Wigner function, Eq. (|l7|), of the relative component ipi P (y r ,t) of the wave packet, Eq. (|lq) , at times i = (a) and 
t — 4-7T (b). The parameters were chosen to be p = 6.0 and 7 = 0.5. 
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FIG. 2: Short-time behavior of an initially Gaussian wave 
packet, Eq. (|l5j), with p = 6.0 and 7 = 0.5. In the upper 
graph we plot the entropy S, Eq. quantifying the entan- 
glement between the two particles. In the lower graph we 
show the interaction energy E-mt, Eq. (|l9|). The £-axis ranges 
from to 2ir, the time of one round-trip of two non-interacting 
particles in a harmonic trap. During the time where the two 
particles interact, seen explicitely via the increase of interac- 
tion energy, the entropy changes, performing a "step" towards 
higher entanglement. 



During times of interaction, indicated by a non- 
vanishing interaction energy E[ n t, Eq. (|l9|), the entropy 
changes considerably. As expected, we have an increase 
in the beginning which goes in a steplike manner. This 
is characteristic for the short time behavior of the model. 



Here the main contributions to 'J/^' are still localized and 
the corresponding Wigner function is not yet smeared out 
over phase space. We emphasize that for longer times the 
entanglement can also decrease during the interaction. 
We discuss this in the next section. 

We have seen that the motion of our two particles is 
dynamically entangled in a steplike process. However, 
the entanglement measure 5, Eq. (|J), is a rather abstract 
quantity, which gives us no hint, how this entanglement 
manifests itself physically. In particular, we shall ask 
whether this state shows nonlocal behavior under certain 
circumstances. In the next section we will look for an 
appropriate situation to answer this question. 

IV. NONCLASSICAL CORRELATIONS 
A. Bell-CHSH Inequalities 

In 1964 John Bell gave certain bounds on correlations 
between two systems within the class of local realistic 
theories (LRT) [|[ These Bell inequalities were later 
formulated in different ways, we will use the description 
by Clauser, Home, Shimony and Holt (CHSH) §, g(J. 

Bell's inequality is based on local measurements per- 
formed on two spatially separated systems i — 1,2. 
Within a LRT the crucial assumption is that the ob- 
served quantities are described by observables Ai(di,X) 
that only depend locally on the measurement parameters 
di chosen for system i and the complete state A of the 
total system. As soon as we define di and A the value of 
Ai is determined. If the simplest case is assumed, namely 
that these values can be only +1 or —1, the observable 
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is dichotomic. Correlations between the observations are 
then given by the quantities (A\(d\, \)A 2 (d 2 , A))lrt ob- 
tained by averaging over all possible states A of the cor- 
responding LRT. The Bell-CHSH inequality now limits 
the range of a combination 

C(dx, d[, a 2 , a! 2 ) = 

Ai(a u X)A 2 (a 2 , A) + A 1 {a' l ,\)A 2 {d 2 ,\) 
+Ai (di , X)A 2 (a' 2 , A) - A 1 (ai , \)A 2 (d' 2 , A) (20) 



of four such correlations to 



\(C) 



LRT 



< 2 



(21) 



for any choice of the four parameter vectors di, d\. 

However, if one describes the observed quantities 
within quantum mechanics, that is by hcrmitian opera- 
tors Ai with eigenvalues ±1 and correspondingly replaces 
the states A by vectors in Hilbert space or density oper- 
ators the correlations of the combination 



C(d\, a[, a 2 , a 2 ) = 

Ai{d{)A 2 {d 2 ) + Aiid'jA^dt) 

+A 1 (d 1 )A 2 (a 2 ) - ii(a 1 )i 2 (a 2 ), 



are now limited by 



K^qmI <2V2. 



(22) 



(23) 



The range between 2 and 2\/2 lies outside the scope 
of any LRT and is therefore a clear signature of the non- 
classical behavior of the underlying system. In the next 
paragraph we shall show that within certain time inter- 
vals also our system can exhibit such nonclassical features 
for an appropriate choice of observables. 



B. Dichotomic Operators for Wave Packets 

It is not obvious how to transfer the assumptions of 
the Bell inequality from two-valued discrete systems to 
continuous systems. Several approaches are conceivable. 
Their distinction comes from the various kinds of two- 
valued observables that can be constructed for a con- 
tinuous system J2l], ^2[ [23] ]. We will rely on the shifted 
parity as our measure quantity as proposed in (^3|. Note, 
however, that this is certainly not the only choice. 

The parity operator acting on subsystem i reads 



Pi 



dxi \ -Xi)(xi\. 



(24) 



This operator has eigenvalues ±1. A parametrization of 
this operator is given very naturally by measuring the 
parity in different reference frames, or equivalently, we 
use the displaced parity operators 



Ai(xi,Pi) = D\{x l ,pi)P i D l (xi,Pi) 



(25) 



for both systems i — 1,2, that is di 
unitary displacement operator 



Di(x,p) = exp 



-{{xpi 



PXi 



(x i ,p l ). The 



(26) 



performs a parallel shift mapping the origin (0, 0) of 
phase space to the point (x,p). If we now choose dif- 
ferent parameters (x\,pi), (x' 1 ,p' 1 ), {x 2 ,p 2 ) and (x' 2 ,p' 2 ) 
we can construct the quantity C, Eq. (fj2"[). 



C. Principles of the Measurement 

The physical interpretation of the dichotomic opera- 
tors Ai(xi,pi) and A 2 (x 2 ,p 2 ), Eq. (p5|), allows us to de- 
scribe a procedure to measure these operators for the two 
particles in the harmonic trapping potential. 

The parity is a global property of a quantum state. 
The operator Pi, Eq. (p4|), seems to demand a quite com- 
plicated measurement process. But if we transform to a 
representation with energy eigenstates \n)i in a symmet- 
ric potential for particle i the parity operator reads 

oo 

/oo 
dxi ^2 \n)i(n\-Xi)(xi\ 

-oo «=° 
oo 

p oo 

= j dx t ^2\n) l {n\x l ){-l) n {x, l \ 

-oo «=° 

OO 

= ^(-lrir^R (27) 



From this we see clearly, that the parity operator is 
dichotomic. A measurement of the parity now means 
to perform a reduced energy measurement that decides 
whether n is even or odd. For a state \ip)i of particle i 
the probability for even parity is given by 



^ (+) =^U(2n|^ 



n=0 



and that for odd parity reads 



P [ t ' =^U(2n + l|V) l 



n=0 



which finally results in 



(Pi) = P, (+) - Pt ] 



(28) 



(29) 



(30) 



The expectation value of the operator Ai(xi,Pi), 
Eq. (|25|), with respect to a state \i/j)i can then be de- 
termined by measuring the probabilities, Eqs. ( |28| ) and 
(p9|), for the displaced state £>i(xi, pi)\ip)i. In addition 
this expectation value i (ip\Ai(xi, pi) \ip)i is up to some 
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FIG. 3: Entropy S (dotted line, right axis) and CHSH sum 
C m ax (solid line, left axis) versus time t for parameters p = 6.0 
and 7 = 0.5. The steps of the entropy, Eq. (|^), already seen 
in Fig. |2| for short times add up to a oscillating long-time be- 
havior. The same can be seen for the sum C max , Eq. jj^), of 
correlations. The bounds for C ma x are given by the value 2 
for any local realistic description and 2\/2 for quantum me- 
chanics. There is clear evidence against a local realistic inter- 
pretation of our model within specific time intervals. 




FIG. 4: Overlap, Eq. @, between an EPR-state \Y r ,P c ), 
Eq. (B4), and the state \^/^(t = 97r)} with parameters p = 
6.0 and 7 = 0.5. The time is chosen such, that the CHSH-sum 
C m ax, Eq. (B2J), takes on its first maximum (see Fig. H). 



scaling the Wigner function W(xi,Pi) of the state 

Hil- 

Since we want to measure the parity of both particles 
simultaneously, we have to separate them in two distinct 
symmetric traps at measurement time t. This prevents 
any further interaction. Then the trap potentials have to 
be displaced [ p2[ , and the reduced energy then gives us 
finally the correlations functions 

(A 1 (x 1 ,p 1 )A 2 (x 2 ,p 2 )) = mA 1 {x 1 ,p 1 )A 2 (x 2 ,p 2 )\^) 

(31) 

for the two-particle state \^f). Note here that in the de- 
scription of this process we have always assumed that the 
two particles are still distinguishable in a certain degree 
of freedom. Otherwise it would be impossible to catch 
them in the appropriate potential. 



D. Violation of the Bell-CHSH inequality 



As stated before, the sum C of correlations, Eq. (22), 
depends on eight real parameters. We used a numerical 
optimization of these correlations within the full param- 
eter space to get the maximal correlations 



C n 



max (C) 



(32) 



for a given state \E'( p -' , Eq. (|l6|). Note that one could also 
think of varying the initial state, Eq. (|l5|), in order to 
achieve even higher correlations. We believe, however, 
that the chosen initial state is the most reasonable one 
if we assume uncorrelated atoms in the beginning. In 
Fig. |^ we show both, the entropic entanglement measure 
S, Eq. (^), and the maximal correlation C max found for 
our model, depending on scaled time t. 

If we look at the entanglement S we still see the single 
steps discussed before. But now they build up an oscilla- 
tory behavior on a much larger time scale and we clearly 



have time intervals in which atomic collisions decrease 
the entanglement. The same oscillations also appear in 
the CHSH correlations. Moreover, they clearly overcome 
the borders of a local realistic description showing almost 
maximal violations of the Ineq. (pl|). 

This strong violation is of particular interest since our 
dynamically generated state is much diff erent from an 
EPR-state as already discussed in Sec. III. This becomes 
clear when we look in Fig. [| at the overlap, Eq. (p"8|), at 
time t = 97r, where we find the first maximum of C max 
(see Fig. |). 



V. CONCLUSIONS 

The prospects of quantum information processing us- 
ing continuous entangled systems clearly demonstrate the 
need to control the generation of such an entanglement. 
Presently, not many simple and controllable systems have 
been described that allow us to start from reasonable ini- 
tial conditions and to arrive at true nonclassical correla- 
tions. We have studied a very elementary model which 
nevertheless leads to an interesting and dynamically gen- 
erated entangled state. This state violates a Bell-type 
inequality when we choose the right set of observables. 
Moreover, it is much different from the ubiquitous EPR 
states which have no simple dynamical equivalent and 
can therefore only be approximated by a physical pro- 
cess with finite resources. In the next step of our study 
we have to understand how these continuously entangled 
states behave under reasonable couplings to an environ- 
ment and we shall devise their potential applications in 
quantum information processing. 
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APPENDIX A: SOLVING THE SCHRODINGER 
EQUATION 

The time independent Schrodinger equation, Eq. (|^), 
can be divided into two parts. If we use the separation 
ansatz ty(y c ,y r ) = ipciy^iprdlr) the equation splits into 
two independent ones: 



dyl 



2E r 



(Al) 



and 



dy 



with constants E = E r + E r 



MVr) = (A2) 



Eq. (Al) for the center-of-mass motion is a harmonic 
oscillator, not influenced by the particle interaction at 
all. The well-known eigenvalues E c = n + i, n <E IN , and 
eigenfunctions 

MEclVc) = (V^2 n niy 1/2 H n (y c ) exp (~y 2 c /2) (A3) 
form a complete basis for this part of the solution. 



Contrarily, Eq. (A2) contains the interaction part 

ffi„t=7%r). (A4) 

The solution therefore splits into two parts 
xp r -(E r ;y r ) for y r < 0, ip r , + (E r ; y r ) for y r > — of 
the equation 



dyl 



Vr 



2E r 



4>r,±{E r ;y r ) = 



together with two boundary conditions 

i> T -(E r ;y r = 0) = tp r ,+(E r ;y r = 0) 

and 



(A5) 



(A6) 



< + {E r ; y r = 0) - $ (E r ; y r = 0) = 27 V» P (E r ; y r = 0) , 

(A7) 

where if) 1 denotes the first derivative with respect to y r . 

Using the ansatz tp r ,±{E r ', y r ) — w{z) exp(— z/2) with 
z = y^ unveils the differential equation 



d 2 w(z) 
''^z^~ 



1 

z 

2 



dw(z) 
dz 



1 - 2E r 



w(z) = 0. (A8) 



w{z) = U I — ,-,z 



(A9) 



which yields 



ipr,±{E r ;yr) = N±U y 4 2gr , ^y 2 j exp (-y 2 r /2) 

(A10) 

with normalization constants A/± . 

To fulfill the constraints, Eqs. (A6) and (A7), we have 
to look for the derivative 



< i± (£ r ;y P = 0)=T2A4 
and the value itself 

ip r ,±(E r ;y r = 0) = Af± 



r((l - 2E r )/A) 



r((3 - 2E r )/4) 



(AH) 



(A12) 



at the origin of space. Then the conditions can be fulfilled 
either by 



N+ = -N- and 



1 



or 



N+ = N- and 



r((3-2£ r )/4) 
27 



(A13) 



T((3 - 2£ r )/4) r((l - 2E r )/A) 

(A14) 

yielding the odd a nd ev en eigenfunctions with eigenvalues 
E r . While Eq. ( |aT^) establishes the undisturbed odd 
eigenvalues 



r ~ 2' 2' 2 '"' 



(A15) 



and eigenfunctions 

ME r ; y r ) = KddH^ r _ 1/2 (y r ) exp (-y 2 r /2) (A16) 
of the harmonic oscillator with normalization 



1 



(i<i 



v 0F2^- 1 /2( jBr -l/2)!' 
the solutions of the transcendental equation 

7 _ 2 

r((3 - 2E r )/4) r((l - 2E r )/4) 

deliver eigenvalues and even eigenfunctions 



/ 1 - 2E r 1 2 

1pr{E r ;yr) = A/ eV en U I ,-,y r 



(A17) 



(A18) 



x exp(-y2/2) 



(A19) 
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FIG. 5: Eigenvalues E r , Eq. (A15) and solutions of Eq. (A18) 



for the Schrodinger equation, Eq. (A2), of the relative motion 
of the two particles with 5-like interation. The interaction 
parameter 7 can be either positive, meaning repulsive intera- 
tion, or negative, corresponding to an attractive force. While 
the odd numbered eigenvalues remain constant with varying 
7, the other eigenvalues increase with stronger repulsion. 




with 



a/;. 



E r 
2 



^(§-%)-V>(i-%) 



(A20) 



where ip(x) is the logarithmic derivative of T(x) 

In Fig. H we show the dependence of the eigenvalues 
of these solutions on the interation strength 7. For no 
interaction (7 = 0) the values correspond to the equally 
spaced eigenvalues of the harmonic oscillator. The en- 
ergy eigenvalues increase due to repulsion or decrease for 
attractive interaction. 

The Wigner function of the ground state is shown in 
Fig. One can see a clear dent in the position distribu- 
tion located at the origin, which indicates the repulsive 
interaction. 



APPENDIX B: EPR-STATES IN POSITION 
REPRESENTATION 

The EPR-State @] for two subsystems 1 and 2 

\X,P) = , jdxi jdx2 exp(iPxi/h) 

2nh, J J 



'2nh , 

x 6{x 1 -x 2 -X)\x 1 )\x 2 ) (Bl) 
fulfills the eigenvalue equations 

[x 1 -x 2 ]\X,P)=X\X,P) (B2) 

and 

\p 1 +p 2 ]\X,P)=P\X,P). (B3) 



FIG. 6: Wigner function of t he g round state of the radial 
wave function ip r (E r ;y r ), Eq. (A19), for an interaction of the 
form Hi n t = 7<5(j/r) with parameter 7 = .5. The energy of 
this ground state, determined by Eq. (A18), is given by E r w 
0.7335, well above the value 0.5 for the non-interacting case. 
We find a perturbed Gaussian distribution having negative 
values along the axis y r = 0. 



If we use Eqs. (|3|) and (|J) we can define scaled EPR-states 

\Y r ,P c ) = -j= Jdy r Jdy c exp(iP c y c ) 

x S(y r -Y r )\y r )\y c ) (B4) 

which are now given in the relative coordinate y r and the 
center-of-mass coordinate y c with parameters Y r = n r X 
and P c = %P. 

We can now compare the state of our two atoms to 
such an EPR-state. We take the state ^^(y c , y r , t) from 
Eq. (|l6|), and find the product 

(y P ,P c |*W(t)) = ^ r (Y r ,t)--j^ Jl Jdy c e^y{-iP cVc ) 



(> M» ( ~2 yc ) 6XP 1^2^ 



ip r (Y r ,t) ' ^lyjexp I --P t 



(B5) 



This leads to the overlap 



(Y r ,P c \¥^(t)) 2 = \MY r ,t)\ 2 • ^exp(-P c 2 ) (B6) 

\ 7T 



9 



which is a Gaussian in the variable P c and reflects the tive coordinate, 
probability distribution of the wave function of the rela- 
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